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Abstract 

A generalized inverse scattering method has been developed for arbi- 
trary n dimensional Lax equations. Subsequently, the method has been 
used to obtain N soliton solutions of a vector higher order nonlinear 
Schrodinger equation, proposed by us. It has been shown that under suit- 
able reduction, vector higher order nonlinear Schrodinger equation reduces 
to higher order nonlinear Schrodinger equation. The infinite number of 
conserved quantities have been obtained by solving a set of coupled Riccati 
equations. A gauge equivalence is shown between the vector higher order 
nonlinear Schrodinger equation and the generalized Landau Lifshitz equa- 
tion and the Lax pair for the latter equation has also been constructed in 
terms of the spin field, establishing direct integrability of the spin system. 



1 Introduction 

Integrable models, since its inception, occupy a distinguised position because 
of their striking predictability and the existence of a special class of extended 
solutions, called solitons. Solitons are characterised by non-dispersive localised 
wave pulses which produce as a result of the competetion between linear disper- 
sion and nonlinearity. Moreover, solitons retain thier shapes even after collisions 
among themselves and consequently exhibit particle like behaviour. Integrable 
models as well as solitons have applications in such diverse areas of physics as 
high energy physics, condensed matter physics, plasma physics, nonlinear optics 
and nuclear physics. 

At present, a handful of nonlinear dynamical equations exist exhibiting soli- 
ton solutions. Historically, KdV equation is the most important one. It 
describes the dynamics of waves moving in shallow water. KdV and KP hier- 
archies also play an important role in describing non-critical strings through 
the construction of r functions. In KdV or KP hierarchy approach, evalua- 
tion of partition function and all correlation functions may be made directly 
through a set of integrable differential equations, bypassing the complicated 
technicalities involved in other approaches. Sine-Gordon equation || is another 
wclknown example of integrable equation, which first arose as a master equation 
for pseudo-spherical surfaces. Besides, it describes the propagation of magnetic 
flux in a Joseption junction transmission line, the propagation of dislocations 
in a crystal lattice and many other physical problems. Nonlinear Schrodinger 
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equation and its higher order generalizations [4,5] have also many applications 
in plasma physics and in nonlinear optics. In particular recently, a lot of ex- 
citement is veered around the soliton solutions of a higher order generalization 
of nonlinear Schrodinger equation, because of its potential application in high 
speed fibre communication systems. The key observation is that the dynamics 
of the higher order nonlinear Schrodinger equation (HNLS) fi not only takes 
care of dispersion loss, but also takes care of the propagation loss as the op- 
tical pulse propagates along the fibre. This is due to the fact that stimulated 
Raman effect Q], which compensates the propagation loss, already exists with- 
ing the spectrum of HNLS equation. As a consequence, a very short pulse can 
be propagetd over a long distance without distortion. Thus, it is optical soli- 
tons, which are responsible behind the successful propagation of optical pulses 
through the fibre over a long distance. Therefore, the importance of the study 
of integrable models and to find their soliton solutions is unquestionable. 

There are several methods exist to obtain soliton solutions of a nonlinear evo- 
lution equation, like Hirota bilinear method Q , Painleve analysis || , Backlund 
transformation Jl(| and inverse scattering transform (1ST) [jn]. 1ST method, 
however, is the most elegant tool, which eventually proves the complete inte- 
grability of the evolution equation. 1ST method, on the other hand, is the most 
complicated one and is intimately connected with the existence of an auxiliary 
linear problem, known as Lax equations. Most of the known integrable evolution 
equations are, in fact, associated with 2x2 and 3x3 dimensional Lax opera- 
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tors. But there are some evolution equations whose associated linear equations, 
in general, cannot be cast in terms of 2 x 2 or 3 x 3 Lax operators. One such 
example is higher order nonlinear Schrodinger equation (HNLS) Q: 

id t q + [3!d xx q + fJ 2 \q\ 2 qie[l3 3 d xxx q + fJ 4 d x (\q\ 2 )q + fJ 5 \q\ 2 d x q] = (1) 

where, q is a complex field and (3i i = 1, 2 • • • , 5 are the constant coefficients. 
Notice that in absense of last three terms, (|l|) reduces to nonlinear Schrodinger 
equation. It is observed recently, for arbitrary values of the coefficients, (0) is 



associated with n x n Lax operators 1 12 . Interestingly, the parameters, (3 are 



related to the dimension of the Lax operators. Unfortunately, 1ST method is 
not well understood for arbitrary dimensional Lax oparators. Nonetheless, 1ST 
methods developed in the context of 2 x 2 [[ll| and 3x3 jl3[ [llj matrix Lax 
operators cannot be generalized straight forwardly to handle more general cases. 
Thus, in order to obtain soliton solutions of ([!]) by 1ST method, a generalized 1ST 
method, associated with n dimensional Lax operator, is to be formulated. In this 
paper, we will, therefore, develop a generalized 1ST method, which ultimately 
leads to solving a set of coupled Gelfand Levitan Marchenko equations jl5| 
for obtaining soliton solutions. We will use the generalized 1ST method to 
obtain soliton solutions for an evolution equation, whose dynamical field is an 
(n — 1) component vector, q. q, in fact, is a vector generalization of the field 
q given in ([!]). The evolution equation may, therefore, be called vector higher 
order nonlinear Schrodinger equation (VHNLS) and will be described in the 
next section. We will also show that under suitable reduction VHNLS equation 



reduces to (|l|). 

The rigidity in the structures of the solitons reveals that the system has a 
huge underlying symmetry. This symmetry is, in general, manifested by the 
existence of an infinite number of conserved quantities. Existence of infinite 
number of conserved quantities is a key criterion for a field theoretical model 
to be called integrable in Liouville sense Construction of Lax pair of a 

dynamical system although itself is a good hint for integrability, many integrable 
models may not satisfy Liouville integrability criterion. In this paper we will 
obtain Riccati equations associated with an n dimensional Lax oparator and 
indentify the generating function of the conserved charges. Consequently, we 
will find conserved charges explicitly in terms of the field variables and their 
derivatives. 

Another interesting aspect, we will address, is to establish a connection be- 
tween nonlinear field theories and spin systems. In this context, it is welknown 
that nonlinear Schrodinger equation is gauge related to Landau Lifshitz equa- 
tion |l7|l . Recently Sasa Satsuma equation |jq] , a particaluar version of HNLS 
equation is shown to be gauge related to a generalized Landau Lifshitz equation, 



where the spin field is associated with ST/ (3) group JL9, 20 . It is thus expected 
that a spin system may also exist corresponding to a vector generalisation of 
HNLS equation. We will establish a connection between the VHNLS equation, 
proposed by us and a generalized Landau Lifshitz equation and obtain a Lax 
pair for the spin system. 
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The organization of the paper is as follows. In section 2, the VHNLS equation 
is introduced and a Lax pair for the nonlinear equation is constructed. We 
also show a reduction procedure under which VHNLS equation reduces to a 
two parameter family of HNLS equation. In section 3, we show the existence 
of infinite number of conserved quantities, by solving a set of coupled Riccati 
equations. We develop the inverse scattering method for an n dimensional 
Lax operators in section 4 and obtain n Gelfand Levitan Marchenko equations. 
Section 5 deals with the solutions ofGelfand Levitan Marchenko equations and 
cosequently find N soliton solutions. We show in section 6 the gauge equivalence 
of VHNLS equation and the generalised Ladau Lifsitz equation and obtain a 
Lax pair for the generalized Landau Lifshitz equation in terms of the spin field. 
Section 7 is the concluding one. 

2 Evolution Equation and Lax Pair 

We propose a vector nonlinear evolution equation of the form 

qt + eqxxx + 3e(g* ■ q x )q + 3e|g| 2 <fx = (2) 

where, the dynamical field variable, q = (g l5 q%, ■ ■ ■ , q n ~i), is an (n — l)-tuple vec- 
tor. The equation (||) is an example of vector higher order nonlinear Schrodinger 
equation (VHNLS). The vector field q may be interpreted as an (n— 1) interact- 
ing optical modes describing the dynamics of a charge field with (n — 1) colours. 
In fact, we will show under suitable reduction, (0) yields to the HNLS equation 




5 



In order to obtain Lax pair for , we first introduce the n xn matrix linear 
eigenvalue problem as 

d x V = U(x,t,A)tf (3a) 

dt® = V(a;,i,A)tf (3b) 

where, ^(x, t) is an n-tuple vector auxilliary field, A is the spectral parameter 
and the Lax operators \J(x,t) and \(x,t) are n x n matrices. Let us now 
assume an explicit form of the Lax pair, XJ(x,t) and ~V(x,t), associated with 
the VHNLS equation, as 

U = -iXS + A (4a) 
V = —eA xx + e(A x A — AA X ) + 2eA 3 (4b) 
-2ieAS(A 2 - A x ) + 4eA 2 A - 4ieA 3 S, 

In the equation (Q) S is a c-no. diagonal matrix and the matrix A(x, t) consists 
of dynamical fields, q{x,t) and q*{x 1 t) only. It is interesting to note that the 
evolution equation for the matrix A(x, t) immediately follows from the compat- 
ibility condition, namely 

V t (x, t) - Y x {x, t) + [Ufa;, i), V(x, t)] = 

provided A and S satisfy the conditions that 

S 2 = l, £A + A£ = (5) 

and as a consequence, the nonlinear evolution equation for A(x, t) becomes 

At + eA xxx - 3e(A 2 A x + A, A 2 ) = (6) 
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It is now clear that various representations of the matrix A in terms of the 
fields q and q* yield to different nonlinear evolution equations. To associate (||) 
with the VHNLS equation (||) , let us consider the explicit expressions of X and 
A(x, t) of the form satisfying the properties (||) as 



n-l 



X — ^ &ii &nn (^ a ) 

z=l 

n— 1 n — 1 

A(x,t) = }^ qj(x, t)e in - }^ q*(x, t)e ni (7b) 

i=l i=l 

where, is an n x n matrix whose only (ij)th element is unity, the rest 
elements being zero and qi(x,t) is the ith. component of the dynamical field q. 

Substituting (^) in (||), the evolution equation becomes 

n — 1 n—1 
qu + tqixxx + 3e(J^ q*j-qjx)q% + 3e(J^ 'ij'lj = (8) 

3=0 3=0 

which is nothing but VHNLS equation (|J), written in the component form. Now 
in order to obtain HNLS equation ([j]) let us consider the following reduction. 
Notice that all the dynamical fields qi in (||) are independent. However, instead 
of (n — 1) independent dynamical fields, if we restrict ourselves to only one 
dynamical field q and its complex conjugate q* in (|^), then all the qiS are not 
independent, g^s, in this case, may be chosen as either q or q* . If we, for 
example, choose m number of qi as q and the rest (n — m — 1) numbers as q* , 
(||) reduces to HNLS equation 

qt + tqxxx + 6me\q\ 2 q x + 3(n - m - l)e{\q\ 2 ) x q = 0. (9) 

The equation (ph, to be precise, is a gauge equivalent version of HNLS equation. 



Soliton solutions of (g) and those of HNLS equation are, in fact, related through 
a U(l) gauge tranformation p"2| . Two well known euqations immediately follow 
from (g). If we choose m = n — 1, (g) gives rise to Hirota equation plj | 

Qt + tqxxx + Ge\q\ 2 q x = (10) 

after rescalling of q and q* as q = (n — and = (n — respectively. 

On the other hand, if we choose m = (n — l)/2, which is only possible for odd 
dimensional Lax pair, (||) reduces to Sasa Satsuma euqation |l8|, |2(], S^j, 

qt + (q X xx + 6t\q\ 2 qx + 3e(\q\ 2 ) x q = 0, (11) 

once again with appropriate scalling of the fields q and q* respectively as q = 
('^yh and q* = (2=i)-$g*. 

3 Riccati Equation and Conserved Charges 

In order to obtain Riccati equation, we first write the Lax equation (^|a,^|aj^) in 
the component form. For the first (n — 1) components of , (ga) can be written 
in the form, 

v lx = -a*! + gi *„ 

^n-lx = -«A* n _i + q n -l^n 
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i.e. 

= -iX^i + qi^n (12a) 

where, i = 1, 2, • • • , n — 1 and denotes the ith. component of But for the 
nth. component of 'J/, the Lax equation (^a) has a different form as 

^nx = «A*„ - ql^i - q^2 ql_{$ n -x 

n-l 

= iWn-^ljVj ( 12b ) 

3=1 

Following now a similar procedure as in ]20| we write, 

r, = |i, (13) 

i = 1, 2, • • • , n — 1, which are related to the conserved charges a nn (X) in the 
following way, 

lna„ n (A) = ln^' n — iXx\ x ^oo 

/oo ^~1 
9? r <) ( 14 ) 

We will see in section 5 that a nn (X) is, indeed, (nn)th element of the scattering 



data matrix and more so it does not evolve with time. By using (12a,b) and 
13[), we may obtain a first order differential equation for each I\, 



n-l 



T ix + 2iXYi - Qj^i ~ 1i = (15) 



3=1 



The set of (n — 1) coupled nonlinear differential equations for 1^ in ( |l5| ) are 
called Riccati equations. It is obvious from (p~4j) that the solutions of Riccati 
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equations eventually determine the conserved quantities. Now in order to solve 
(p^|), we assume a series solution of r, as 

oo 

r i (x,A) = ^Q(x)A-" (16) 

n=0 

Substituting ([IB]) into (p5|), the following recursion relations may be obtained. 

CI = 0; C* = | (17a) 

and 

fc-l n-l 

2*2+2 + ((7l +1 )x - E C fe-i+m E = ( 17b ) 

m— J — 1 

with k = 0, 1, 2, The infinite number of Hamiltonians (conserved quan- 

tities) may explicitly be determined in terms of the dynamical field variables qi 
and their derivatives by expanding a nn (X) in the form, 

Zna„„(A) = (n-l)E7^2lTT^ A ^ ( 18 ) 
and thus comparing ([l8]) with ([TJ) and (|l6|), Hi becomes 

hi = (_!)/(„ -i) y ^ ^ 

The explicit expresseions of the first few low order Hamiltonians are given below. 



1 f 

Hi = — - / dx[E g^'] (20a) 

1 /" ™ _1 

^ 2 = ^ — rj / ' Ar X '/>/" ( 20b ) 

- „ n— 1 n— 1 

#3 = — T / Gb[(E ffi«J*) 2 - E ( 20C ) 

« - 1 ™ 1 



j=l J'=l 
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^ „ n—1 n—1 

Hi = , _ jr / dx^2[q*q jxxx - q* xxx + 3(^ \qk\) 2 {q*jqjx - q* jx qj)} 

^ ' J j=l k=l 

(20d) 

1 » n—1 n — 1 n—1 

h 5 = ^ T /^[E^- + 2(E^%) 3 -(E(i^i 2 )-) 2 

n—1 n—1 n — 1 n—1 

fe=i j=i j=i fe=i 

We verify directly by using equations of motions that Hi, H2, H3, Hi and H^ 
are, indeed, constants of motions. 

Let us now specialize to HNLS equation. If we assume that m number of g^s 
are chosen as q and the rest n — m — 1 as q* , the Hamiltonians in ( pp| ) reduce 
to 

Hi = J dx\q\ 2 (21a) 
# 2 = (2m - n + 1) / dar(?*g* - gg*) (21b) 
Hz = j dx((n - l)\q\ 4 - q x q*) (21c) 
#4 = (2m-n + l)J dx[3(n-l)\q\ 2 (q x q* -qq*) + ^(q xxx q* -qq* xxx )} 

(21d) 

H 5 = J dx[q xx q* xx + 2(n-l) 2 \q\ 6 + [(2 m -n-l) 2 -4(n-l)}\q\ 2 q x q* x 

[(^-l)+ 2m( ( ^:^ 1) ]((M 2 )x) 2 ] (21e) 



Notice that -Hi in (|2l|a) has a universal form, which implies that all solitons 



11 



have the same energy irrespective of their shapes. However, the momentum, 
H2 crucially depends on the numbers of q, chosen in a given representation 
of the matrix A ([7]). For example, H2 becomes zero for Sasa Satsuma case, 
i.e. for m = (n — l)/2. In fact, all conserved quantities having even indices 
become trivial for Sasa Satsuma case. H3 in (|2l|c) also has a somewhat universal 
form, depending only on the dimensions of the Lax pairs. But the higher order 
conserved quantities starting from H5 do not possess such universal forms. Their 
explicit forms depend both on the dimensionality of the matrix Lax pair and 
also on the representations of the matrix A. 

4 Generalized Gelfand Levitan Marchenko Equa- 
tions 

We now generalize 1ST method suitable for studying n dsmensional Lax opera- 
tors. The first step in this direction is to obtain a set of generalized Gelfand Lev- 
itan Marchenko equations. This generalisation is a nontrivial one and crucially 
depends on the properties of scattering data matrix. However, we will see that 
for three dimcnisonal Lax oparators generalized Gelafand Levitan Marchenko 
equations will reduce to Sasa Satsuma case jl8| as is expected. We broadly 
follow the treatment of Manakov, developed in the context of 3 x 3 Lax oper- 
ator |l|] and for that assume Jost functions, for real A, satisfy the boundary 
conditions, 

$« = e ie - lXx (22a) 
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with i = 1, 2 n — 1, but the nth. one satisfies a different boundary condition 

like 

$(«) = e n e iXx (22b) 

as i — > — oo. Similarly, as i — > oo, other set of Jost functions satify the 
boundary conditions 

¥« = ei e- iXx (23a) 
for i = 1, Ti — l and for the nth. one, 

#0) = e„e iA:E . (23b) 

In the equations ( |22| ) and ( p3| ) e^s are the basis vectors for an n-dimensional 
vector space. It follows from (^b) that = —A for real valued A and thus we 
have 

d x (*(i) f *(2)) =0 (24) 

for any pair of solutions of equation (@a), and VlK 2 ), having the same 

eigenvalue. It is straightforward to show from dzj) and ( f23| ) that 

$Wt$0) = $Wt$W) = ^ (25) 

for i,j = 1, 2, n. Since the set of Jost functions are linearly independent 

and the maximum number of independent Jost functions is n, we may express 
the set of Jost functions $, as a linear combination of 'J, as 



$( i )(i,A)=5;ay(A)#(i I A) (26) 
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where ay (A) is the (ij)th element of scattering data matrix, which can be 
expressed by using (|25|) and (^6|) in the form 

ay (A) = ^(x,X)^(x,X) (27) 

The orthogonality property of the scattering data matrix elements for real eigen- 
value A, subsequently follows from (^2|), (|2|) and ( p7| ) and thus we obtain 

n 

2^ ctifc (A)ajfe(A) = 5ij (28) 
fc=i 

which finally gives 

n 

A) A) (29) 

j'=i 

It is further interesting to see, by exploiting the properties of [ay] and [ay] 
matrices, that we can write the element a*, as the cofactor of the elements of 
the matrix [ay]. In particular, a* ni element can be written as 

<i = (-l) n+i det[a ni ] (30) 

where [a n i] is a (n — 1) x (n — 1) matrix, constructed from the n x n scattering 
matrix, [ay] with nth row and ith column being omitted, i.e. det[a n i] is the 
minor of a n i element of scattering matrix, [ay]. Now by using ( p6[ ) and (|30|), 
we obtain the following useful relations among the Jost functions $W and 
The first n — 1 Jost functions in ( p6|) satisfy 



-i n- 1 * 
_ ^ J2( Ad ^n))kj^ j) e iXx = *W e iXx - _2*$<») e <** (31a) 

with k = 1, Ti—l, but the nth. Jost function, $w obey the relation 

1 1 

$W e -^ = ^Wg-iA* + ^ ohj-tf^e - *** (31b) 
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Notice that in deriving (|3l|), we have used the following properties of the 
scattering data matrix. 

n-i 

a *nk&ij = ^ [«nfck {Add [<Xnk])lj 

1=1 

It is important to mention that analyticity properties of the Jost functions and 
the elements of scattering matrix may be obtained from (|3l|). 

We are now going to derive the Gelfend, Levitan and Marchencho equation 
for an n dimensional Lax pair. Let us consider an integral representation of the 

Jost function vljW for (i = 1,2, , ri). For the first n—1 Jost functions we may 

choose the following integral representations 

/>oo 

[x, A) = e je - iXx + / dyK u) (x, y) e -* Xv (32a) 

J X 

with j = 1, 2, n — 1, while the nth Jost function may be written as 

*(")(a;,A) = e n e lXx + / dyK^ (x, y)e lXy . (32b) 



where, et, i = 1, 2, • • • , n are basis vectors for an n dimensional vector space 
and the kernels K^-* and K*-™) are n dimensional column vectors, which may be 
written explicitly in the component form as 

n 

K^(x,y)=J2 K( dH^V>m (33a) 

m— 1 

n 

KM(x,y) = Y,K£ ) (x,y)e™ (33b) 

rn—1 
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Substituting ( p2| ) in (31b), we obtain 

_L$W = e n e iXx + / di/KW^i/Je^ + ^^-e^e-^ 

+ E — / ^ (j) (^)«" iAv (34) 
Multiplying now both sides of (34) by i f^° oo d\e~ lXz , with an assumption 
z > x and using the analyticity property of the Jost function, it follows that 

2 7 rl 00 dA a„„(A) e _ M + A M 2tt ^ a„„(A) ^ 
+ rPTdy g ^ K« (x, „) e -^0 (35a) 

The L.H.S. of (||a) can be simplified further by taking into account that t-tt 

is analytic in the lower half plane except at the points, say A* with ImX* > 

and j — 1,2 N, where - 1 ^ has N simple poles. Moreover, we assume 

that at the simple poles A*, to be of the form 

n-l 

^(x,\*) = J2c^^ p) (x,X*) (35b) 



With these assumptions, L.H.S. of (35 a) becomes 



3 - =1 ^nnK^jJ p=1 

where a^ n denotes derivative with respect to A. By Substituting the integral 
representations of *W(x,Aj) from ( pa) ), ( |35^ ), i.e. the L.H.S. of (§^a) finally 



reduces to 

N p-iA*z "-1 poo 

- l T,7^1T^T, C nJ^ X ' jX + dyK^(x,y)e^y] (35d) 
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Let us now introduce a funtion F p (x + y) as 

F(x | y)=i Y C ^^ +V) , f (36) 

In terms of the function, + y) in (|3^), (|35|a) and ( |35d| ) together may be 

written in a compact form as 

n— 1 w— 1 „oo 

K(")(a;,z) + ^e p F p (x + z) + ^ / (x, y)F p (z + y) = (37a) 

p=i p=i •' x 



The inetgral equation (37a) is one of the desired Gelfand Levitan Marchcnko 
equations for the kernel K^™). Other integral equations for the kernels K^) for 



p — 1, 2, • • • , n — 1 may be obtained from (31a) and (R2) in a similar way like 



that of (37a). The integral equations for the kernels K.^ thus turn out to be 
of the form 



K^{x,z)-e n F;{x + z)- / dyK {n \x,z)F;(z + y) = 



(37b) 



provided z > x. In deriving (37b) we have used the following identity 



Cn P = = ^[anp{\])]ki{Adj[a np {\j)]) a 8 k i- 



The set of coupled equations (p7f) may be called as generalized Gelfand Levitan 



Marchenko equations. Substituting now (33b) and (|37b|) in (37a), to a first 
approximation, we find the Gelfand Levitan Marchenko equation for the pth. 
component of Kl"': 

n—l />oc />oo 

KW(x, z) + F p (x + z)+Y, / dsK^(x, s) / dyF m (y + z)F* m (y + s) = 

rn=l Jx Jx 

(38) 

which will be used later to find the soliton solutions for VHNLS equataion. 
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5 N Soliton Solutions 



To obtain soliton solutions, let us associate dynamical fields qi(x,t) with the 



kernels, if. in ( |38[ ). Substituting (|32b| ) into Lax equation (|2|a) we find 



qi (x) = -2K\ n \x,x) 



(39) 



and consequently, it is evident that the solution of ( |38[ ) gives rise to soliton 
solutions in terms of scattering data elements. But before going to solve (|38|), 
we first compute time evolution of scattering data element. Notice that as 
\x\ — > oo, the Lax equation (3b^b) leads to 



~dt 



-4ieA d £ 



(40) 



which, in turn, determines time evolution of the scattering data elements. For 
example, scattering data elements, a nj) (A, t) for j — 1, 2, • • • , n — 1 evolve with 
time as 



a nj (X,t) = a n j(\,Q)e' 



while the element a nn (X,t) is time invariant: 



(41a) 



a„Ti(A, t) = a„„(A, 0), 



(41b) 



which eventually justifies our conjecture in section 3 that a lm , indeed, can be 



associated with the conserved quantities. From (|29[) and ( 35b ) it follows that 



the coefficients C„p also satisfy the similar time dependence as a n j in J41a|) and 



thus 



C^(A,t) = C«(A,0)e- 8 ^ 3t 



Tip 



(42) 
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If we now restrict ourselves to soliton sector, a np (A) becomes trivial and the 
function F p (x + y), in this case, reduces to 

j— i fin V l j ) 

Time dependence of the function F p (x + y) may be obtained immediately from 
© as 

^(* + tf)=*2^ — — Vta^ ■ (43) 

To solve the integral equation ( |38| ) for the soliton solutions, the kernels 
K p n \x + y) are assumed to be of the form 

N 

K^(x + y) = J2^Me- iX ^ (44) 

3=1 



Substituting (43) and (44) in (p8|), we obtain a set of n — 1 algebraic equations: 

N n—1 —i\*x 

K«>(x,z) + F p (x + x)- J2 E 



Cnr (0)Cnr fe (0) ;( 2Afc _ A ;_ Anx+8fe (A|-Af )t = n ^ 

<„(A*R'„(A fc ) 1 ] 

Solving ( |45| ) for K p (x 1 x) and subsequently using (^), the TV soliton solutions 
for each dynamical field qi(x,t) may be expressed as 

N 

94 (M) = -2E( BCrl )y e ~ lA> ( 46 ) 

where, B and C are respectively (n — 1) x N and N x N matrices whose explicit 
forms are given by 

(Bk^C^We- 8 ^ 3 '-^ 
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Let us now consider an explicit form of one soliton. One soliton solution, 



which follow from (46), may be written in the following form 

*(M)= '^f l X l 2 i (47) 

1 + Ep=l \ F p( X + X )\ (Ax-A-V 

where 



qi{x,t) in ( |47| ) may be expressed in a more conventional way, by choosing the 
position of the pole at Ai = |(— £ + zr/) and by introducing e 7i+Jl5i 
and thus 



2r ) < Il (A*) 



e 



7 I -» ) a;+er,( J 7 2 -3C 2 )t+i(<5 i +^+e4 3 t-3e4 2 r ) t) 



* } = 4 ^ l + ^-^e^W^)* (48) 



which can further be simplified as 



by introducing 



qi(x,t) = —t — 49 

coshP 



e 2r„ = y- e 2 7p 
p =i 

Pfoi) = V x - e V { V 2 - 3£, 2 )t - T„ 
Q<(a;,i) = ^ + <(£ 2 - 3£r/)< + <f< 



Once again, if we specialize to HNLS equation, each independent field qt{x, t), 
depending on the models, reduces either to q or to q* and as a consequence qi 
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in (El) yields to 



= /, ' hA (50) 
V(n — ljcos/m 



where, i = r/x - e^rj 2 - 3£ 2 )t - 7 - ±/n(n - 1) and B = £x + e£(£ 2 - 3£rj)t + S. 
It is interesting to note that we have obtained precisely the same expression for 
one soliton in |fl2| . 

6 Generalized Landau Lifshitz type equation as 
the Gauge equivalence system 

We now show an interesting connection between the VHNLS equation and the 
generalized Landau Lifshitz type equation by exploiting the gauge equivalence of 
the Lax pairs of these two dynamical systems. The procedure is similar to that 
between the nonlinear Schrodinger equation and the standard Landau Lifshitz 
equation ]l7j . 

Under a local gauge transformation, the Jost function, ^(x,t, A) changes to 

* = g- 1 (x,t)*(x,t,\) (51) 

where g(x,t) = ^(x, t, A)|a=o- We claim that g(x,t) is an element of SU(n) 
group. As a consequence of the gauge transformation (|5l|), the Lax equations 
(U|) become 

* x = U(x,t,A)* (52a) 
* t = V(x,t,A)* (52b) 
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where U and Vave the new gauge transformed Lax pair, given by 

V(x,t,X) = g- 1 (V-V )g (53a) 

V(x,i,A) = g- 1 (V-V )g (53b) 

with U =U| A=0 = g x (x,t)g~ l (x,t) and V = V| A=0 = g t (x,t)g~ 1 (x,t). This 
leads to 

A = g x {x,t)g- 1 (54) 

Since A belongs to su(n) algebra, g{x,t) obviously belongs to SU(n) group, 
which justifies our claim. We may now identify the spin field of the Landau 
Lifshitz type equation as 

S = g- l {x,t)Y.g{x,t), S 2 = l (55) 

With this identification, the gauge transformed Lax pair ( pT| ) may be expressed 
in terms of the spin field S ( |55| ) and its derivatives only, yielding 

U = -iXS (56a) 
V = -4ieX 3 S + 2e\ 2 SS x +ie\{S xx + ^SS 2 x ) (56b) 
In deriving ( p6|) we have used the following important identities 
SS X = 2g- 1 Ag 
SS 2 X = -Ag-^Aig 
S xx + SS 2 x = 2g- 1 VA x g 
The zero curvature condition of (|56|), namely 

u t - V x + [U, V] = 
22 



ultimately leads to the generalized Landau Lifshitz type equation 

3 3 

St + eS xxx + —e(S x + SS XX S X + SS X S XX ) — (57) 
withS G SU(n)/(U(n-l)) n . 

7 Conclusion 

We have formulated a generalized 1ST method for an n dimensional Lax opara- 
tor. Subsequently, the generalized 1ST method is used to obtain N soliton 
solutions for a multi-component generalization of HNLS equation, proposed by 
us. We have, although, obtained soliton solutions for VHNLS equation, 1ST 
method developed here is quite general and is applicable for obtaining soliton 
solutions for all integrable nonlinear equations. The sech structure of one soli- 
ton solution of VHNLS equation is quite interesting, particularly in the context 
of nonlinear optics, since it can easily be produced from the output of a mode 
locked laser. 

We have shown the integrability of VHNLS equation also in the Liouvillc 
sense. This has been acheived first by finding a set of coupled Riccati equa- 
tions and subsequently by identifying the generating function for the conserved 
charges. It is found that the last diagonal element of the scattering data ma- 
trix may be identified as the generating function of the conserved charges. 
This is also confirmed from the time evolution of scattering data matrix ele- 
ments. We have also established an intriguing relationship between the VHNLS 
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equation and a generalized Landau Lifshitz equation, where the spin field S £ 
SU (n) /(U(n — 1))". Moreover, we have obtained a Lax pair for the spin system 
implying direct integrability of the generalised Landau Lifsitz equation. 
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